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Abstract 

The angular power spectra of cosmic microwave background are an- 
alyzed under the light of the evolutional scenario of the universe based 
on the renormalizable quantum theory of gravity in four dimensions. The 
equation of evolution is solved numerically fixing the power law spectrum 
predicted by the conformal gravity for the initial condition. The equation 
requires to introduce a dynamical energy scale about lO^^GeV, where the 
inflationary space-time evolution makes a transition to the big-bang of 
the conventional Friedmann universe. The quality of fit to the three-year 
data of WMAP implies the possibility to understand the observation by 
quantum gravity. 
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Observations of the cosmic microwave background (CMB) anisotropics by 
the cosmic background explorer (COBE) [I] and the Wilkinson microwave anisotropy 
probe (WMAP) [51 (3] Sj have established the inflationary scenario of the uni- 
verse. Various cosmological parameters have been determined precisely based 
on the cosmological perturbation theory [71 Hj ■ In those analyses, however, 
the primordial spectrum close to that of the Harrison-Zel'dovich [5| is assumed 
as the initial condition of the big bang scenario of the universe. One of the 
important problems we wish to discuss in this letter is the dynamics to produce 
such a primordial spectrum in terms of the fundamental point of view rather 
than introducing an ambiguous field from phenomenological purpose. 

It has been recognized that before the big bang there is an inflationary 
expanding epoch which solves the horizon and the flatness problems [TOl [TTj . 
The inflation naively says that the universe grew about the order of 10^" from 
the birth to the present, which implies that the size of the Hubble distance 
today was placed within the Planck length at the Planck time. Thus, it is 
natural to consider that the primordial spectrum was originated from quantum 
fluctuations of space-time and these Planck scale phenomena are recorded as 
the CMB anisotropics observed by COBE and WMAP [12]. 

Since the observed spectrum of constant l{l + 1)C/ suggests an almost scale 
invariant primordial spectrum for the large angle correlation, it is natural to 
expect that initial fluctuations are generated by a dimensionless field, and the 
only candidate is the metric field within known physical fields in the microscopic 
scale. 

The inflationary scenario induced by quantum effects of gravity was first 
proposed by Starobinsky in 1979 [llj. His idea has an advantage that it is not 
necessary to introduce any additional fields. At that time, however, the idea 
cannot explain how the inflation terminates to the big bang after expansion with 
sufficient number of e-foldings and why tensor fluctuations are less dominated 
than scalar ones in the very early universe while quantum gravity effects are so 
strong. 

The renomalizable quantum conformal gravity we employ here can explain 
these problems in terms of the asymptotically free property of the traceless 
tensor mode in the metric field whose dynamics is governed by the Weyl action. 
It suggests that there is a dynamical energy scale separating between quantum 
and classical space-time, and beyond this energy scale, an inflationary quantum 
space-time is realized, and then the conformal- mode fluctuation dominates. In 
this way, Starobinsky's idea of inflation is revived on the foundation of the 
modern field theory. 

Quantum fluctuations of the conformal mode are getting small during the 
inflation. Applying the linear perturbation theory about the inflationary back- 
ground, we compute the transfer functions from the Planck time to the big bang 
which took place at the dynamical scale [T3]. We then calculate the multipole 
components of the CMB angular power spectrum comparing with the WMAP 
data. 
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Evolutional Scenario of The Universe Renormalizable quantum gravity 
[14l[T5] is defined by tiie four-derivative conformal invariant actions, {—^/t^)\/—gC^ 
and —b^/^G4, where C^^act is the Weyl tensor and G4 is the Euler density, in 
addition to the Einstein-Hilbert action and conformal invariant matter fields. 
The dimensionless coupling constant t in the Weyl action is introduced to take 
care of the traceless tensor mode at the short distance scale, while the constant 
6, which is introduced to renormalize divergences proportional to the Euler den- 
sity, is not an independent coupling constant because the Euler term does not 
have the dynamical components. Quantization is carried out perturbatively 
about conformal flat space-time with the vanishing Weyl tensor, and thus the 
metric field is expanded as gf^i, — G^'^g^u with g^^ = rj^^ + h^iy + • • where </> 
and are the conformal mode and the traceless tensor mode, respectively. 

The renomalized coupling, tr, for the traceless tensor mode is shown to be 
asymptotic free, whose beta function was computed in [HI [14] as /3 = —(iot^ 
with /3o > 0. This justifies the perturbative treatment for this mode, and also 
implies the existence of a dynamical scale Aqg, where the running coupling 
constant is written as 

lA2(p)=/3oln(pVA|G) (1) 

for a physical momentum p. The asymptotic freedom yields that the coupling 
constant is getting small at very high energies, and configurations with the 
vanishing Weyl tensor are chosen quantum mechanically that prohibit a singular 
configuration with a divergent Riemann curvature tensor. 

The conformal mode is quantized non-perturbatively so that the conformal 
invariance becomes exact when the traceless-mode coupling constant vanishes 
at very high energies. The dynamics of the conformal mode is induced from the 
measure as the Wess-Zumino action of conformal anomaly, known as the Riegert 
action, containing the kinetic term (— &i/87r^) x y/^cfjA^cj) [T71[TSl[Tni[2n], where 
A4 denotes the fourth-order conformal invariant operator defined on g^^. The 
coefficient has been computed within the lowest order as hi = (27Vx -I- lliVw + 
1247VA)/720-f 769/180 19J, where Nx, A^v and Na are the numbers of conformal 
scalar fields, Weyl fermions and gauge fields, respectively. 

Evolution of the early universe is divided into three stages by two mass scales, 
the reduced Planck scale Mp = I/vSttG and the dynamical scale Aqg, ordered 
as Mp 3> Aqg [H]. At very high energies beyond the Planck scale, the space- 
time is dominated by the quantum fluctuation with exact conformal invariance. 
The symmetry begins to be broken about the Planck scale, toward the stage of 
inflationary expanding universe with the Hubble constant Hp, = y^Sn'^/biAIp. 
The running coupling gradually increases during inflation, and it diverges at 
the dynamical scale Aqg- It is the period that the conformal invariance is 
completely broken and correlation length becomes short to result the emergence 
of classical space-time. 

The dynamics of the inflationary phase before the big bang is effectively 
described by including corrections of the traceless-mode coupling constant to 
the Wess-Zumino action such as &i(l — oii^ -I- • • •) = biBQ{tr) with oi > 0. 
The higher order effects are taking into account by a resummation form Bq = 
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1/(1 + ait^/n)'^, where k is a parameter that hes in the range < k < 1. 
Under the spirit of the mean-field approximation, we simphfy the momentum 
dependence of the running coupling by its time-dependent average: replac- 
ing the physical momentum by the inverse of the proper time r as l/t^(T) = 
/3o ln(l/r^AQQ). It shows that the running coupling diverges at the dynamical 
time scale 1/Aqg (= ta), and then the dynamical factor Bq vanishes indicating 
the transition from the conformal gravity to the Einstein gravity. In this way, 
we obtain the evolutional homogeneous equations of motion [131 115j , 

Bo(r) [h +7HH + AH^ + l^'^H + ei/"*) - 3iJg [h + 2H'^^ ^ (2) 

and the conservation equation 

Bo(t) {2HH -H^ + 6H^H + SH^'j - 3HIH^ + Sn^p/bi = 0, (3) 

where H is the Hubble parameter defined hy H ^ a/a and a = e"^ is the scale 
factor. The dot denotes the derivative with respect to the proper time and p is 
the matter density. 

The space-time initially evolves in an infiationary expansion with H = Hj^. 
The number of e- foldings from the Planck time I/Hd (= rp) to ta is approxi- 
mately given by the ratio of two mass scales: Afe = log[a(TA)/a(rp)] ^ Hu/Aqg, 
which will be set about Hd/Aqq = 60 as popularly accepted. The coefficient 
of the Wess-Zumino action bi is taken as 15 and 20 for the analyses to compare 
to the WMAP3 data in the following. The dynamical scale is then given as 
Aqg ^ lO^^GeV. 

The other parameters in the model are rather insignificant. Since they de- 
pend on the non-perturbative dynamics of the traceless mode, they are chosen 
phenomenologically as (3o/bi = 0.06, ai/bi — 0.01, and k — 0.5. With these 
combinations of parameters the homogeneous equation ([2]) and the evolution 
equations for scalar and tensor fiuctuations, which will appear later, preserve 
their forms independent of bi. The number of e- foldings is then computed to be 
Afe = 65.0. The solution of the evolution equations is depicted in figure [T] for 
the case of &i — 15, where Hj^ is normalized to be unity. The sharp increase of 
the matter density at the transition point indicates the big bang where energies 
stored in the conformal mode shift to the matter degrees of freedom. 

Below the energy scale Aqq, the Einstein action becomes dominant and the 
space-time makes transition to the classical phase. This phase is described by 
the low energy effective theory of gravity expanding in derivatives of the metric 
field |13| in an analogy to the chiral perturbation theory for QCD. For simplicity, 
we connect the conformal universe to the Friedmann universe at the transition 
point. This simplification will not cause much effects in later discussions because 
we consider fiuctuations with the size of the Planck length at the Planck time. 
Therefore, at the transition point the size is extended much more than the 
correlation length, and the patterns of spectra are insignificant on the dynamics 
and the parameters Po, ai, k employed at the transition. 
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Primordial Spectra Since the inflationary solution is stable, gravitational 
fluctuations about this solution get smaller during the inflation. Let us first 
give a rough estimation for the amplitude of scalar fluctuation considering a 
dimensionless contrast of the scalar curvature fluctuation, 5R/R. The denomi- 
nator is the curvature of inflationary background, which is 127?^ for de Sitter 
curvature with H = Hj^ . Since the curvature has two derivatives, the curvature 
fluctuation would be order of the square of the energy scale, E"^. Hence, the 
amplitude of scalar fluctuation is estimated to be SR/R ~ E'^/12H'^. This im- 
plies that we can apply the linear perturbation theory for the density fluctuation 
under the inflationary background in the period from the Planck time rp to the 
dynamical time t\. At the dynamical energy scale, the amplitude of the scalar 
fluctuation is estimated as AQQ/12ffQ ~ l/12J\f^, which gives the magnitude 
about the order of 10~^ similar to the observations. 

For the scalar fluctuation we compute the evolutions of the so-called Bardeen 
potentials defined by ds^ = a2[-(l -I- 2'^)drf + [l + 2<^)dy?]. The coupled 
equations for the evolutions of these fields in the inflationary background have 
been derived in [T3]. The initial configuration at the Planck time is given by 
$ = 4" because of the dominance of conformal-mode fluctuations at very high 
energies, while at the transition point the dynamics requires the configuration 
satisfying $ = — The scalar spectrum defined by the two-point correlation 
of the Bardeen potential is given by the quantum conformal gravity as 

ppi(/c) = A, (fc/m) (4) 

at the Planck time, where A: is a spatial comoving momentum and m = a(Tp)ifD 
is the comoving Planck scale at the Planck time. Since this scale appears in the 
evolutional equations of fluctuations, it is a dynamical parameter to determine 
the pattern of the spectra at the transition point. The dimensionless amplitude 
V^Aj is given by the order of 10~^ obtained by substituting E ~ H-o into the 
expression of curvature fluctuation estimated above. The scalar spectral index 
is given by the anomalous dimensions of the scalar curvature [211 HI] , 

= 5 - 8(1 - Vl-2/6i)/(l - (5) 

In the large hi limit, it approaches to the Harrison-Zel'dovich spectrum with 
Us = 1-h 2/61 + 4/6? -fo(l/6?). 

The tensor spectrum at the Planck time is given by 

Pf{k)^At(k/mYK (6) 

Because of the asymptotic freedom for the traceless tensor mode, the amplitude 
is considered to be much smaller than the scalar amplitude {At <C A^) and the 
index is given by = 0. 

The transfer functions in the inflationary period are defined by $(rA,A:) = 
^(ta, T'p)<i>(rp, /c) for the scalar mode and /i^'^(ta, fc) = Tt{T\, Tp)hJj^{Tp,k) for 
the tensor mode. The primordial spectra at big bang are obtained as Ps{k) = 
Z^{TA,Tp)PP\k) and Ft(fc) =r,2(rA,rp)pP'(fc). 
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The transfer functions are computed numerically, and it is shown that the 
amplitude of the scalar fluctuation gets smaller as estimated above, while the 
tensor fluctuation stays to be small [13]. The scalar spectrum that is initially 
blue (ris > 1) shifts to red in higher momentum regions with k > m. From this 
behavior we can fix the comoving Planck scale. Since the WMAP data favors 
the scalar spectral index less than one at fe = 0.05Mpc~^, the value of m is 
taken to be as small as 0.04Mpc~^. In figure[2]we show the result of simulation 
for the case of bi = 15 and 20 with m = 0.04Mpc~^. This value of m implies 
that the scale factor at the Planck time is a(Tp) ~ 10^^^ when we consider the 
size of Planck length at the Planck time grows up to the size 1/m today. 

CMB Multipoles The CMB angular power spectra are calculated using the 
cmbfast code [35]. The initial conditions to be supplied are chosen from the 
transfer functions obtained at the transition time specified above. The coeffi- 
cient of the Wess-Zumino action bi which depends on the matter contents is 
taken to be larger than that for typical GUT models0 so that the scalar spec- 
tral index is close to one for the momentum region less than m. In figures 
[3] and m we show the numerical results of TT, TE, EE, and BB spectra for 
bi = 15(n,, = 1.15) and 20(71^ = 1.11) together with the WMAP3 data. The 
other primordial parameters such as the tensor-to-scalar ratio r and the ampli- 
tudes are adjusted to fit the observation data. In order to see the dependence 
of the primordial parameters, we employ the cosmological parameters to be the 
same as the best fit values of the analysis of WMAP [3j . 

The suppression of the low multipole components may be explained as the 
consequence of the dynamical scale. If we wish to discuss it rigorously we need to 
evaluate the two-point correlation function for long-distance separation, which 
will involve non-pcrturbative dynamics of the traceless mode. Instead, we here 
simply give a phenomenological ansatz based on the dimensional analysis. The 
damping factor for a spatial separation far beyond the correlation length is 
assumed to have the following form: k"^ /{k'^ -\- uX^), where A = a(rp)AQG = 
m/60 is the comoving dynamical scale at the Planck time. This factor should 
be multiplied to the initial spectra ^ and (O. The simulation results with the 
damping factor are depicted by the dashed lines in figure [3] and IH There is 
small improvement on the I — 2 component in the TT spectrum, but it seems 
to suggest that we need further investigations on the dynamics of the two-point 
correlation requiring more sharp fall off. 

For higher momentum region more than k > 2m, we need to study with 
considering non-linear effects of the evolution equation, because the fluctuation 
with high momentum beyond this range corresponds to that with the momentum 
far beyond the Planck scale initially. For the scalar fluctuation, the non-linearity 
of conformal field theory becomes significant until the amplitude reduces small 
enough to apply the linear approximation safely. We regard the non-linear effect 
hopefully sustains the amplitude so as not to fall off at high momentum region 
which is required by other observational results. As for the tensor fluctuations, 

iPor 5(7(5), SO(IO), and Ee models, bi ~ 9, 12, and 18, respectively 
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the linear approximation is applicable in the high momentum region due to the 
asymptotically free property of this mode. 



Conclusion The renormalizable quantum gravity based on the conformal 
gravity in four dimensions suggests that there appears a dynamical energy scale 
Aqg about the order of lO^^GeV, between the Planck scale and the GUT scale. 
This energy scale separates the conventional classical space-time from quan- 
tum space-time with conformal invariance. The big bang can be understood 
as a space-time transition where energies stored in extra degrees of freedom in 
higher-derivative gravitational fields shift to matter degrees of freedom. 

The CMB angular power spectra arc calculated on the basis of the quantum 
gravity. We numerically simulate evolutions of the scalar and tensor fluctuations 
during the inflationary period from the Planck time to the transition time, and 
by setting the computed spectra at the transition as the initial conditions of the 
cmbfast code we obtain the TT, TE, EE, and BB angular power spectra. The 
results fit to the WMAP3 data, which suggests that the renormalizable quantum 
gravity is one of the possibility to explain the origin of CMB anisotropics. 

The condition that the dynamical energy scale is being lower than the Planck 
mass implies that quantum effects turn on much larger scale than the Planck 
length, and thus very high energy particle is dressed by quantum gravity and 
space-time is deformed locally to prevent the particle itself from forming a black 
hole. Such a quantum gravity effect about the dynamical scale might be ob- 
served in future experiments such as observations of gamma-ray bursts [23' , in 
addition to CMB. 
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Figure 1: The time evolution of H and p for bi = 15 and Hy) /Aqq = 60, where 
Hj) is normalized to be unity. 



8 



bi=l5,m=0M 
Z?i=20,m=0.04 



0.5 




' ^ 

0.05 0.1 

momentum k [Mpc"^] 

Z7i=15,m=0.04 

Z7i=20,m=0.04 




momentum k [Mpc' ] 



Figure 2: Patterns of the scalar and tensor spectra, Ps and Pt, for bi = 15 (thick) 
and 20 (thin). In each case, the simulation is carried out by m = 0.04Mpc~^, 
and the amplitude is normalized appropriately. For the tensor spectrum, there 
is no significant dependence on the value of bi. 
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Figure 3: The TT speetra for bi = 15(thick) and 20(thin) with m = 0.04Mpc"^ 
The tensor-to-scalar ratio r and the damping factor u are determined to fit 
with the WMAP3 data. The amphtude is normahzed at the first peak. The 
part at low-multipole components is enlarged and depicted to see the r and u 
dependence for each 61, while there is no significant dependence on u at high- 
multipole components. The other cosmological parameters are fixed to be the 
best fit values in [3] such as Qh = 0.041, f^cdm = 0.205, f^vac = 0.754, Hq = 73.1, 
Te = 0.108, Tcmb = 2.726, and Fhc = 0.24. 
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Figure 4: The TE, EE, and BB polarization spectra. 
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